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ABSTRACT 


This thesis investigates the multiplexing of Mach-Zehnder type 3x3 terminated fiber¬ 
optic sensors demodulated by either 'quadrature' or 'symmetric' methods using intensity 
modulation of the source. 3x3-couplers produce signals that permit unmodulated passive 
demodulation of interferometric signals. The theory is described and the results of a 2x1 
element array optical demonstration are presented. Possible architectures using this 
demultiplexing technique are presented for several applications with different return line 
requirements. The technique was successfully demonstrated and warrants further 
investigation to increase the number of sensors and reduce the number of return lines for 
specific applications. The multiplexing technique presents the opportunity for possible cost 
savings over other phase generated carrier techniques, which require wavelength 
modulation of the source, and significant optical path differences in the interferometers, 
and are therefore constrained to presently very expensive sources. The technique 
presented uses compatible low coherent laser sources such as Compact Disc quality (830 
nm) devices. 
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TABLE OF SYMBOLS 


Symbol 

Definition 

A 

The sum of the squares of the two interfering electric fields and will be termed 
the dc bias value of the optical intensity. 

B 

Twice the product of the two interfering electric fields which reprents the 
amplitude of the cosinusoidally varying signal whose argument is the 
differential phase.optical intensity in the fiber, i.e., the maximum variation from 
the dc value. 

C 

The amplitude of the mixing sinusoidal signal. 

c„ 

The coupler ratio for coupler n The coupler ratio is the power split ratio of a 
particular coupler as a fraction of the total input power, i.e., a percentage or 
equivalent fraction. 

F 

The number of return lines for a given array architecture. 

G 

Gain factor. 

Jn 

The Bessel function of order n. 

k 

The wave number in a vacuum defined as k = 2nlX, where X. is the 
wavelength in meters unless stated otherwise. 

L 

The physical length of optical fiber in meters unless stated otherwise. 

Lc 

The coherence length of a laser, defined as the distance over which the 
amplitude modulation falls off to an arbitrarily determined fraction of the 
complete modulation value [Reference 6, page 46]. 

M 

The number of lasers in an array of sensors. 

m 

The modulation index, a measure of the amount of modulation as a fraction of 
the total available, having a value between 0 and 1. 
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Symbol 


Definition 


//, V, and w 


The number of sensors per laser in an array of sensors. 


The refractive index of a particular medium where the refractive index is 
defined as the ratio of the speed of light in a vacuum to the speed of light in the 
medium. 


Pf, denotes the output power from a particular laser. Subscripts are used to 
denote the power from a particular coupler in an array, e g., the power 
transmitted by the coupler. P' denotes the power normalized to Pq. 


The differential phase shift in the arms of the interferometer. In general ^ is a 
time dependent quantity, however the time dependent notation is dropped 
for compactness. It is often used with various subscripts to denote the origin of 
the phase shift e g., <l>(jrij^t) denotes the phase shift due to factors other than 
the signal of interest. denotes the amplitude of the sinusoidal signal of 
interest. 


The amplitude of wavelength modulation which occurs at the modulation 
frequency. 


The total number of sensors for a given array. 


The frequency of the phase shift in the interferometer, either in Hertz or 
Radians respectively. 


The light intensity of each of the three legs of a 3x3 terminated Mach-Zehnder 
interferometer. Each is phase shifted 120° relative to the other two. 
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I. INTRODUCTION 


.1. MOTIVATION 

The United States Navy is seeking to take foil advantage of fiber optic technology, in 
particular, the bringing together of fiber optic communications technology and fiber optic 
sensor technology. Fiber optic communications systems have many advantages over other 
methods of data communication. These include high bandwidth, low power consumption, 
security, electromagnetic interference (EMI) protection, electromagnetic pulse (EMP) 
protection, light weight, low loss, and low cost. Optical fiber sensors based on 
interferometry have been designed with high sensitivity and dynamic range Various 
demodulation and multiplexing schemes are feasible for use with these interferometric 
sensors [References 1, 2, 3, 4, and 5]. For example it has been demonstrated that 3x3- 
coupler terminated interferometric sensors with quadrature [Reference 6] and 'symmetric' 
demodulation yield particular advantages [References 7, 8, and 9]. An array of these 
optical fiber sensors in an optical architecture has the potential to accommodate both the 
benefits of the interferometric sensors and optical communications systems. An all-optical 
system removes the need for outboard electronics and thus has the potential for great 
weight and cost savings. Such an optical system using 3x3 couplers and 'symmetric' 
demodulation is demonstrated in this thesis. 

B. GOAL 

This thesis demonstrates that the use of 3x3-coupler terminated interferometric 
sensors and either quadrature or 'symmetric' demodulation enables the use of amplitude 
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modulation (AM) frequency-division-multiplexing (FDM) techniques where the laser 
intensity is the quantity that is amplitude modulated 

The theory and performance of such a multiplexing scheme is investigated with sensor 
demonstrations using both analog electronics and optical fiber interferometric sensors On 
the basis of this research a design for an array of M x N sensors is presented which may 
provide large cost savings with improved system performance. 


C BACKGROUND 
1. Interferometric Sensors 

Interferometry is a technique used to measure differential phase shifts The optical 
phase delay generated in a physical length of fiber L is. 

= hiL, ( 11 ) 


where n is the refractive index of the fiber core and k is the wave number in a vacuum (k - 
2x/X). Neglecting birefringence from the two orthogonal polarizations supported by the 
fiber, the fractional variation in phase can be expressed as; 


d<t> 


k 


dn 


n 


dL 
L ’ 


( 12 ) 


where the changes in n and L are related to a particular measurand (e g., pressure, 
temperature, magnetic field, or other environmental changes). 

Various techniques have been demonstrated to selectively enhance the change in 
phase in the fiber for a particular measurand (e g., acoustic waves), while reducing 
changes in phase due to unwanted effects (eg., acceleration). Mach-Zehnder 
interferometric sensors use two-beam interferometry as illustrated in Figure 1.1. It has 
been shown that by using both legs as signal legs, in push-pull configurations, improved 
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sensitivity and, through design, highly favorable properties (e g., acceleration canceling ) 
can be obtained. 



Figure 1.1. Schematic diagram of a Fiber Optic Mach-Zehnder 
Interferometric Sensor 


An ideal 3x3-coupler provides three signals of equal amplitude, each with a relative phase 
difference of 120°. The intensities of each output leg of the 3 x 3-coupIer are: 

u = A 5cos + (13) 


Bcos + + 

(1.4) 

Bcos 

(1.5) 


where A is the value of optical intensity from the fibers (analogous to a dc offset), B is the 
peak value of the intensity and is the static phase shift resulting from the difference 

in optical path length between the two signal legs. This quasi-static (or slowly varying 
compared to the signal of interest) phase also includes the phase shift due to 
environmental effects (hydrostatic pressure, temperature, etc.) as well as the physical 
difference in the length of the legs. Examples of these types of signals are shown in 
Appendix A. The signal of interest <j)(t) is assumed to be of the form; 
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^(0 = ^,sin(<y,0 


(16) 


Using the trigonometric identity: 

cos(a+>S) = co^(a)cos(^+sin(a)sinO!?), (17) 

the intensity, equation (1.3), can be rewritten as; 

M = ^ + 5cos[^,sin(a>,/)]cos[^^^(/)+^,] 

- fisin [<!>, sin {a, /)] sin [^rf„y,(/)+ 

Similar expressions for v and w with phase shifts of ±120°, as appropriate, also can be 
written. This equation is of the form of a Bessel function generating function. The 
generating function and associated series from Reference 10 (page 361, equations 9 1 42 
and 9.1.43) are reproduced below; 

oo 

cos(zsin ff) = Jq (r) + 2^*/„ (’) cos(2A:^ (1.9) 

k=\ 


sin (r sin 0) = ,(z)sin[(2A: + l)4 0 10 ) 

k=0 


The intensity, equation (6), can therefore be written as a series of Bessel weighted sine and 
cosine functions as follows. 


u = A + Bcosi(^j„^it))\ 




lr=l 


( 1 . 11 ) 


-^sin(^,,„^(/)) 


2^ (^. )sin [(2A: + l)m,/] 

. *=o 


Thus the frequency spectrum consists of components at multiples of co^ which are in turn 
weighted by Bessel function terms and the drift term, Bcos{^j„j,(t)). Appendix A plots 
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the lower order Bessel functions and the frequency spectrum of Equation 1.11 for typical 

HO and H„fi(0 

2. Symmetric Demodulation 

The technique of symmetrically demodulating the outputs was developed by 
Professor Robert M. Keolian of Naval Postgraduate School and is detailed in References 
7, 8, and 9. A schematic of the scheme is shown in Figure 1.2. Modifications and 
performance specifications of the symmetric analog demodulator of Reference 9 are 
detailed in Appendix B. 



Figure 1.2. Passive Symmetric Demodulation Scheme (adapted from Reference 8). 










The demodulation scheme is based upon the following expression; 


4>(t) 


□C f u(v - w) + v(w - u) w(u ~ v) 
^ (u^+v^+w^) 


(1 12 ) 


where ^i) is the signal of interest to be recovered and the dots above a variable indicate a 
time derivative. The implementation detailed in Reference 9 is outlined below for 
continuity. 

The output intensity of each leg of the interferometer is as previously defined; 

u = A + Bcos[^K0 + ^</„/j(0]. (113) 

with similar expressions for v and w. After photodetection and amplification by low-noise 
transimpedance amplifiers, the symmetric voltages are added together and scaled by a 
factor of -1/3 to produce -A. This average dc content signal is then summed to each of 
the symmetric, interferometric voltage signals to remove the average dc component, A. 
These interferometric signals, with the average dc component removed are denoted again 
by II, V, and w, where 

u = 5cos[^/)-i-^^„^(/)] 

V = ficos[<^(/)+^,„/t) + 120°] (1.14) 

w = 5cos[$K0 + ^,fr,y» (0-120°]. 

These three symmetric interferometric signals are cross multiplied by the difference 
of the 'complimentary' derivatives to yield; 
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u{v-w) = 5cos(^)|^5^sin(^-120‘’) + 5^sin(^ + 120°)j 

v{w-u) = 5cos(^-120° sin(^ + 120®) +5^ sin(^)j (115) 

w(m-v) = ficos(^-120“)|^5^sin(^ + 120“) + 5^sin(^)j. 

It can be shown by trigonometric identities or with phasor diagrams that the first equation 
of (1.16) can be simplified to; 

m(v-h') = y/3B^ ^ cos^(^). (1 16) 

By noting: 

cos‘(x) + cos*(jf + 120*) + cos‘(x-120°) = 4, (1 17) 

the sum of the three equations in (1.16) can be shown to result in: 

Ji5'[cos'(^) + cos'(^+120^) + cos2(^-I20®)] = (1.18) 

The factor of B^ varies as a function of laser intensity, temperature, and the polarization 
angle of the light in the fiber. This dependence is removed by summing the squares of each 
of the interferometric signals (with the dc component removed), resulting in: 

(1.19) 

and using this result as a normalizing division factor, the ratio of Equation 1.18 to 
Equation 1.19 is then integrated to produce 

Unity scale factors have been assumed on all the operators (summers, dividers, etc.) 
in presenting the demodulation technique. 
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Following on from the background of interferometry and demodulation, the theory 
and demonstration of multiplexing interferometric sensors using intensity modulation of 
the source is introduced. 
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li. MULTIPLEXING OF INTERFEROMETRIC SENSORS 


The proposed multiplexing scheme uses 'Optical Fiber Subcarrier Multiplexing' 
(OSCUM) [Reference 1. pages 544-5, and Reference 3, pages 167-8] although this 
approach has not previously been demonstrated with interferometric sensors. Reference 2 
(section 4.5, pages 167-8) details the use of this multiplexing technique with intensity 
sensors. 

A. THEORY 

The theory will be presented using Figure 2.1 and applying subscripts to the equations 
which relate to the numbered 'position' in the Block Diagram. 



Figure 2.1. Multiplexing/Demiiltiplexing Block Diagram: Circled numbers relate 
to subscripts in Equations 2.1 - 2.9. 
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The laser is driven with an AM modulation current given as. 


4 ,(/) = + msin(to.r)]. 

and the resulting light intensity from the laser I,(t) is then: 

A^{t) = i4,[l+ /nsin(tu^/)]. 

The electrical signal from the output of the sensor, at point 2, is then given as 

= [>4 + 5cos(^)][l+/nsin(n;„/)] 

= A + Bcos{^) + mA sin(fi;„/) + mB cos(^)sin( 6 >„/), 

where the notation of Reference 11 has been used. Here (j> denotes our signal of interest 
and includes the drift and coupler terms, given as; 

^ =^O+^d„y,(O + ^c = ^,sin( 6 ),t)+^j„^(/) + /(120°), (2.4) 

where / = 0, +1, or -1. Using the trigonometric identity: 

sin(a)cos(^ = ^sin(a-)3) + -Lsin(a+^), (2 5) 

the signal at point 2 can then be rewritten as. 

v^{t) = A + Bcos{^) + mA sin i(o„{) + -^[sin (o}^ r - ^) + sin (o)^ t+ <!>)], (2.6) 

where the gain as a result of the transimpedance amplifier and fiber losses are all 
encompassed in the amplitudes of each term. An example of a computer generated 
spectrum of the signal represented by Equation 2.6 is shown in Figure 2.2. The peak signal 
phase modulation ^5 = ti radians, the drift phase term is zero, the signal frequency is 1 
kHz and the modulation frequency is 30 kHz, and the intensity modulation depth m = 
0.75. Figure 2.2 clearly shows the baseband and modulated portions of the signal. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 
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Figure 2.2. Frequency spectrum of an interferometric output phase modulated 7C 
radians peak, at 1 kHz, and amplitude modulated at 30 kHz. 

The signal v^O), given by Equation 2.6, is now mixed with a scaled voltage 
proportional to the laser drive current modulating signal. This mixing signal has no dc 
term and is defined as: 

Vo(/) = Csin(<y„/). (2.7) 

Thus the signal at point 3 becomes: 

Vj{t)=AC sin(<y„/) + 5Csin(<y,,/)cos(^) 

V M. O ) 

+ mAC sin'{o)j) + mBC sin^(<y„/)cos(^). 


which can be rewritten as: 


11 















































(29) 


mAC mBC /a 

Vj(0 = —- + cos(<^) 


BC 

+ AC sia((2}j) +-[sin(<y,/ + ^>) + sin(^y^/ - ^)] 

+ ^d£.cos(2(U^/) + — ^ —[cos( 26>^/ - ^) + cos(2<u^/ + ^)], 


where the top and bottom lines of Equation 2.9 represent the modulated portion of 
Equation 2.6 (or Figure 2.2) frequency translated to baseband and twice the carrier 
frequency, respectively. The second line of Equation 2.9 represents the baseband portion 
of Equation 2.6 (or Figure 2.2) translated to the carrier frequency. Using our previous 
example (Figure 2.2) and mixing with a signal of frequency of 30 kHz (C = 1 from 
Equation 2.7), the spectrum shown in Figure 2.3 is obtained. Figure 2.3 clearly shows the 
translation caused by mixing and how subsequent recovery of our signal of interest by low 
pass filtering can be achieved. It will be shown later that this mixing or heterodyning is 
necessary because the modulated signals will overlap in the baseband, precluding their 
recovery by simple filtering. 

The other two legs (relatively phase shifted by 120° from each other) are similarly 
mi.xed and filtered and can then be demodulated using the passive symmetric demodulation 
scheme previously described. 
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Figure 2.3. Spectrum of signal of Figure 2.2 mixed with a signal, phase and 
frequency locked to the laser modulating signal at 30 kHz. 


B. DEMONSTRATION WITH OPTICAL SYSTEM 

It was verified that the ratio of components in the baseband was the same as the 
modulated band (neglecting the DC and carrier components) by use of an analog 
electronic simulation of the optical components and signals. This indicated that the static 
phasing would be preferred with this approach. Appendix C describes this 
demonstration and contains all the significant results obtained. The essential result of the 
demonstration was the verification of the computer generated simulation of Figure 2.2. 
Once this demonstration was successfully completed, a two-laser, two-sensor optical 
demonstration was tested; the schematic is shown in Figure 2.4. Two lasers are amplitude 
modulated at two different frequencies and respectively. The two sensors 
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inleirogated by these lasers are then coupled together before photodetection, 
demultiplexing, filtering and demodulation. 



Figure 2.4. Scheinatic for Optical Demonstration. 

A computer simulation example of two sensor signals and two modulation 
(subcarrier) frequencies is shown in Figures 2.5 and 2.6. They follow directly from the 
derivation of the previous section but highlight the potential for overlap of sideband 






signals and the heterodyning technique used to recover the sigruil of interest This example 
represents the multiplexing and demultiplexing of the interferometric signals from one of 
the three legs shown in Figure 2.4. Figure 2.5 shows two modulated interferometric 
signals with the following characteristics; subcarrier modulation frequency = 40 kHz, 
sensor signal amplitude = l .Stt radians, sensor signal frequency = 1 kHz, bias angle 
(pji = 45°, and f „2 = 60 kHz, (f >,2 = radians,/,^ = 600 Hz, <I>J 2 = 0°. The amplitudes and 
bias angles of the two sensor signals were specifically chosen to be different to highlight 
the overlap at baseband. 



Figure 2.5. Example of two modulated interferometric signals. 

Figure 2.6 shows the signal of Figure 2.5 after heterodyning with modulating signal at 
frequency The overlapped baseband signals have been translated to f„, (40 kHz) 
and the interferometric signals of interest can be recovered by lowpass and bandpass 
filtering. The recovered interferometric signals from each of the three legs can then be 
demodulated. The filtering must provide the same phase response to each of the three legs 
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(ideally constant over the frequency range of interest) so that the interferometric signals 
can be accurately demodulated. 



0 20 40 60 80 100 


Frequency (kHz) 

Figure 2,6. Signal represented in Figure 2.5 heterodyned with modulating signal 
at frequency (40 kHz). 


16 






















































1. Experimental Setup 

The equipment setup for (he optical demonstration is shown in Figure 2.7. 
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The driver and modulation section uses the Seastar LD200 Ultra Stable Power 
Supply to drive the semiconductor lasers. This driver can modulate the intensity up to a 
frequency of about 1 MHz. The optical power output starts to drop off at about 800 kHz 
The variable voltage divider in Figure 2.7 is required because of the conflicting 
requirements of the AD633 multiplier and the external modulation input of the LD200 
The external modulation input of the LD200 converts a voltage input to a laser drive 
current output. The ratio of the driving voltage to resulting drive current is 10 mV to 
1 mA [Reference 12]. Thus a voltage of the order of 50 to 100 mVpp is required to 
provide the required drive current modulation levels. The AD633 multiplier chip, on the 
other hand, prefers to see the largest signal possible (up to 10 V); so, by use of the divider, 
a phase matched signal can be simultaneously provided to the LD200 (50-100 mVpp) and 
theAD633 (10 Vpp). 

Figure 2.8 shows the detection, transimpedance amplifier, and amplification for one 
of the three legs. The two gain stages are included to amplify the signal before mixing so 
that a reasonable level signal (comparable to the laser modulating signal) is seen by the 
multiplier in the mixing stage. The amplifiers remove the DC component of the signal and 
each provide a noninverted gain of 20. The 20 kQ resistor in the amplifier input path is 
required to provide a return to ground for the very small input current [Reference 14, 
section 3.05, page 95]. The noise up to this stage and introduced by the amplification will 
also be amplified, i.e., there will be no increase in signal to noise ratio (SNR). This was not 
considered an issue in this demonstration of the multiplexing technique but would have to 
be addressed in an actual sensor system demonstration. 

Figure 2.9 shows the mixing and low pass filtering circuit where THB292 is an 
active filter integrated into a single package [Reference 13, page 118] with a 3 dB cutoff 
at 22 kHz. The filter has a frequency response as shown at Appendix D. This low cost 
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integrated package ($12 85 each) was chosen to simplify the implementation. In a 
practical system the cutofT required may well be far greater than 22 kHz. A rule of thumb, 
developed by Dr. David A. Brown and extensively verified by experience, gives the 
interferometric signal bandwidth as a function of signal frequency and amplitude as. 

A = 5^./,. (2.10) 

where the subscript sh denotes sigiuficani bandwidth. This estimate is considered 
conservative but, nevertheless, provides a reasonable estimate of the maximum frequency, 
signal level, or combination of both that can be demodulated by the system if we designate 
the cutoff frequency of the low pass filter as X* • 
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Figure 2.8 Detection and Traiisimpedniice Ampliner Stage. 


For example, a system with a dynamic range of 100 dB at 20 kHz with an assumed 
minimum detectable signal of 5 prad/^Hz and a maximum of 4% THD (a typical value 
using Demodulator Version 2.1 [Appendix B]) would tolerate a maximum signal 
calculated as follows; 
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f 

20 log 

V 


/ffad/VH z ^ 
5 /ffad / >/Hz ^ 




100 dB re Imrad / VHz at 20 kHz 
0.5 rad. 


and a lowpass filter cutoff above: 


( 2 . 11 ) 


= 5x0.5x20x10’ = 60kHz. (2.12) 

This cutoff frequency then allows for a maximum signal at 100 Hz, and a maximum of 4% 
THD, of 

^ = A = = 100 radians. (2 13) 

5/, 500 

or a dynamic range of 

= 20x10*^ = 146dBrel/irad/VHzatl00Hz, (2.14) 

5x10-^ 

assuming the same noise floor. 

The benchtop interferometers of Figure 2.7 are described in detail at Reference 9, 
pages 60-66 (see also Appendix D). The purpose of the benchtop interferometer is to 
provide a signal from a push-pull optical sensor that can be easily adjusted by a function 
generator connected to the benchtop interferometer. The three outputs of the circuit of 
Figure 2.9 are connected to the 9207 resistor of the analog synunetric demodulator, 
version 2.1, described in detail in Appendix B. This allows the transimpedance amplifier of 
the demodulator to be reconfigured as a voltage amplifier and saves having to rewire or 
reconstruct the demodulator front end, depending on whether or not the demodulator is 
being used with the demultiplexing scheme or the direct demodulation of interferometric 
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sensor signals. This stage oF amplification is still used to ensure that the three signals have 
the same ac amplitude, just as in the direct demodulation configuration [Appendix B]. 
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Figure 2.9. Circuit of the mixing and low pass filteriug stage of Figure 2.7. 

2. Results 

Prior to coupling the two benchtop interferometers together, each was tested 
separately. A modulation signal was applied to the appropriate laser and then the sensor 
signal was demultiplexed and demodulated. Figures 2.10 and 2.11 show the demultiplexed 
demodulated outputs for inputs to both channels prior to coupling together. 
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200 }.isec/div. 


Figure 2.10. The upper signal represents the input to benchtop interreroineter ^1 
(<|>,i == n radians amplitude) while the lower signal represents the 
demultiplexed/demodulated output (f_, - 80 kHz). 









500 nsec/div. 

Pigtire 2.11. Upper sigiinl represents the input to beiichtop interferometer U1 
(({>,2 = 2.5 n rndi.'ins amplitude) while the lower signal represents the 
dcniultiplexed/deniodulated output (r„, 2 ’= 120 kHz). 









After connecting the two benchtop interferometers, with a 50/50 coupler, the 
simulated acoustic signals for each interferometer were again demultiplexed and 
demodulated. The outputs for two different inputs to each of the benchtop interferometers 
are shown in Figures 2.12 and 2.13. 
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Figure 2.12. Interferometers 1 and 2 coupled together. The upper 
signal represents the input to benchtop interferometer Hi (4»,i = 2.5 n 
radians amplitude) while the lower signal represents the 
demultiplexed/demodulated output (f^, = 80 kHz). 
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Figure 2.14 shows the demultiplexed and demodulated output for a triangular wave 
benchtop interferometer input (2.5 V at -1350 Hz). Even though the frequency spectrum 
of triangle wave has numerous components above its fundamental (1350 Hz) and the low 
pass filter cutoff frequency is 22 kHz, the demodulator is able to demodulate the band- 
limited demultiplexed system. This figure was chosen to highlight the robustness but 
ultimate limitations of the technique. 
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Figure 2.14. Interferometers 1 and 2 coupled together. The upper signal 
represents the input to benchtop interferometer ffl while the lower 
signal represents the demultiplexed/demodulated output (f^, » 80 kHz). 
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Figures 2.IS and 2.16 show the baseband and demultiplexed signals, respectively, 
for a 500 Hz benchtop interferometer signal. The ratios of the components of the signals 
appear to be in the same ratio but Table 2.1 shows they are not. The reader is reminded 
that this was not the case when the multiplexing scheme was tested with analog 
electronics. This difference in the ratios of the components may be due to a number of 
causes (or a combination thereof) These include the phase response of the low-pass filter 
not being constant over the entire passband, the optical path difference (OPD) in the 
interferometer legs and the drift angle change that occurred between plots. Independent of 
the reason for the change in ratios of components, the change is the same in all three legs 
of the interferometer and does not affect the ability of the demodulator to correctly 
demodulate the signal. This would tend to indicate that the cause is due to the OPD in the 
interferometer or the change in drift angle, rather than the phase response of the filters 
Further investigation of the cause of this change was not carried out since it was always 
the same in all three legs and essentially amounted to a drift angle change between the 
baseband and demultiplexed signal. 

Although Figures 2.15 and 2.16 do not show the noise floor, the measured 
difference between noise floors in the baseband spectrum (-114 dB re 1 V/VHz at 5 kHz) 
and the demultiplexed spectrum (-98 dB re 1 V/VHz ) was 16 dB at 5 kHz (see Appendix 
D, Figures D8 and D9). A noise floor of -98 dB re 1 V/VHz equates to minimum 
detectable signal of 13 ttrad/VHz. This is a large difference but is seen to be due to the 
electronics since the optical crosstalk between sensors was less than 60 dB. The crosstalk 
was measured by mixing with f„i and measuring the size of the known spectral 
components due to interferometer #2. The electronics for this demonstration was 
assembled on a bread board with many leads to and from the demodulator. This 
preliminary measurement of the noise floor is considered only an estimate of the system 


27 







noise floor The increase in noise floor and hence detectability threshold due to the 
nuiltiplexing scheme should not be more than about 3-4 dB as a result of the modulation 
losses and demultiplexing electronics. 
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Figure 2.15. Baseband speclriim of 500 Ilz benchtop interreroineter 
signal (Note 60 Ilz noise) prior to multiplexing, = 0.6 radians. 
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Figure 2.16. Spectrum of signal in Figure 2.15 after demultiplexing. 
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Table 2.1 Frequency component ratios for Figures 2.15 and 2.16. 



C. PRACTICAL LIMITATIONS 

As a result of the optical demonstration several important practical limitations of the 
multiplexing scheme became evident. These considerations did not have any significant 
effect on the demonstration of the technique but would be important considerations in a 
complete system demonstration. The most important of these is semiconductor laser diode 
characteristics and the subsequent system bandwidth limitations. 
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1. Semiconductor Laser Diode Characteristics 
A plot of the optical power output ver^s input drive current curve for an ideal laser 
diode would indicate that the maximum amplitude modulation and corresponding bias 
point utilize the complete linear range of laser operation [Figure 2.17], Figure 2.18 shows 
an operating curve that has been exaggerated to highlight the nonlinearities and power 
saturation curvature. 



Figure 2.17. Bias point and amplitude modulation range for a 
laser diode (after Reference 16). 


The 'kinks' or nonlinearities in the operating region are a result of inhomogeneities in 
the active region and power switching between different lateral modes of the laser. The 
downward curvature at the upper end of the operating curve, or power saturation, is 
attributed to active layer heating [Reference 16, page 105]. 
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Diode curreni 

Figure 2.18. Example of kinks and power saturation in laser diode 
operating curve (after Reference 16). 

Details of the causes and methods to reduce these nonlinearities is well documented 
(e g., Reference 1, chapter 5; Reference 5, chapter 3; and Reference 16, chapter 4). 
Operating curves for all the lasers used in the demonstration are included (Appendix D) 
but are far too crude to appreciate the significance of these nonlinearities. These 
nonlinearities cause higher harmonics of the modulation frequency. Modulation of sensor 
information at higher harmonics of the carrier represents lost energy that cannot be 
recovered by the demultiplexing and demodulation process. 

The nonlinearities can be better quantified by measuring the frequency spearum of 
the optical output of an intensity modulated laser and comparing this to the fi-equency 
spectrum of the modulating drive current. Figure 2.19 shows the frequency spectrum of 
the HP3314A function generator used to generate the intensity modulation of the lasers 
shown in Figures 2.20 and 2.21. The spectrum of the HP3314A is within specifications 
[Reference 17, page 74], 
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The two lasers shown are typical of semiconductor laser diodes with the second 
harmonic approximately 20-30 dB down from the fundamental being the typical range 
stated in several of the references. One of the best traditional techtuques to reduce this 
harmonic distortion to 30-40 dB involved the circuit technique of quasi-feedforward 
compensation [Reference 18] but this technique has been superseded with the advent of 
special purpose highly linear devices. The linearity of such devices is measured in terms of 
the third harmonic intercept (TOI) [see Appendix D, Figure D7]. These lasers are more 
expensive ($750-$1200 US) but may be required depending on the complete system 
requirements. It should also be noted that these lasers can only be obtained in the 1300 
and 1550 nm wavelengths. 

Also, if less expensive sources are used, the higher order harmonics of the carrier 
impose a bandwidth limitation on the system. Modulation of different sources has to be 
carried out in one of two ways; either at frequencies between the lowest carrier and its 
second harmonic, or at carriers separated from each other by more than twice the 
frequency of the next lower carrier frequency. The latter method can be used since the 
third and higher harmonics are generally negligible. Both these methods impose a stricter 
limit on the number of possible sensors than with the use of highly linear sources where 
higher order carrier components are negligible. 
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Figure 2.19. Frequency Spectrum orHP3314A Function Generator for 
desired output of 10 Vpp at 100 kllz used to drive lasers of 
Figures 2.20 and 2.21. 
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Figure 2.20. Frequency SpecCrum oniifaclii IfL83l2G #SP27f2 Inser 
driven by Senstnr LD200 power supply fed by signal of 
Figure 2.19. 
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Figure 2.21. Frequency Spectrum of Sharp LT015MDO ttSP2676 laser 
driven by Seastar L0200 power supply fed by signal of 
Figure 2.19. 
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2. Aging and Coherence Length Considerations 
The coherence length of the laser must be greater than the optical path difference in 
the interferometer legs. Using the typical data shown at Reference 19, page 6S. a 
coherence length of approximately 7 mm could be inferred. This was calculated using a 
full line width dA of 0.1 nm (estimated from Reference 19) and an operating wavelength of 
830 nm from; 



This Lq of 7 mm was presumed to be close to the optical path difference of the benchtop 
interferometer. The causes and effects of aging of semiconductor laser diodes are well 
documented [Reference 1, Volume I, section 5.3.3.4, Reference 5, page 56, and Reference 
16, section 4-5]. One cause is related to the drive current level. The laser had been 
operated at high drive levels and this may have possibly reduced the coherence length of 
the laser. The laser was replaced with a Hitachi HL8312G which was presumed to have a 
longer coherence length (i.e., greater than the OPD in the interferometer), if only because 
it was new. The fringe visibility returned. This problem highlighted the importance of 
coherence length and optical path difference considerations and, especially, their 
relationship as a function of aging due to high drive current levels, high temperature 
operation, or any other aging mechanism which will change the operating characteristics 
of the laser diode. 

3. Operating Wavelength Considerations. 

The choice of 830 nm as the operating wavelength for the optical demonstration of 
this thesis was motivated by the inexpensive CD type semiconductor laser sources 
available at this wavelength. It is, however, very important to consider the overall cost of 
the various system components and the direction in which the optical fiber 
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communications industry and technology are moving. Fiber for 1300 nm sources is now 
cheaper than that for 830 nm. The highly linear source described earlier are only available 
in the 1300 and 1550 nm ranges. Optical amplifiers, which may be a consideration in 
multiplexing (considered later), are presently only available in the 1550 nm range. 

All these cost and component availability issues need to be considered when 
choosing the operating wavelength for an all optical sensor array. It may well be that the 
developments of the communications industiy will determine the most economical 
operating wavelength. In either case this multiplexing technique provides a less expensive 
source option. This technique indicates the possibility for use with inexpensive CD type 
sources at 830 nm ($150-350 US), or relatively inexpensive, highly linear sources similar 
to those for cable television applications at 1300 or 1550 nm ($750-1500 US). These 
sources are both inexpensive compared to the highly coherent sources required by other 
presently used PGC FDM techniques involving wavelength modulation of the source [e g.. 
Reference 2, pages 170-172], in which sources typically cost $10,000-15,000 US. This 
cost saving is significant enough to warrant further research into applying the technique to 
larger arrays at RF frequencies. With the present cost of the sources, fiber, demodulator 
and demultiplexing, the most expensive component in the system is the 3x3-coupler in 
each sensor which presently cost $400-750 US each, while 2x2-couplers cost on the order 
of $100 US. 
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in. DESIGN OF AN M X N ARRAY 


A. MULTIPLEXING ARCHUECTURE FOR AN M x N ARRAY 

The multiplexing technique demonstrated optically and described in Section II B, 
utilizing two lasers and two sensors, can be extended to an A/ x ^ sensor element array. 
Figure 3.1 shows a possible implementation topology for a 3x4 element array. Each of 
the bold vertical lines connecting a stave of sensors (or vertical column in Figure 3.1) 
actually represents three separate fibers and has only been drawn as one for clarity in the 
figure. The three interferometric signals obtained for each stave contain the sensor 
information from all the sensors in that stave. To reconstruct the individual sensor 
information the three interferometric lines must be demultiplexed once and low-pass or 
band-pass filtered as appropriate before demodulation. Alternatively each sensor's 
interferometric signals could be demultiplexed and low-pass filtered prior to 
demodulation. For either approach to demultiplexing, a demodulator is required for each 
sensor in the array if ail are to be required to be interrogated simultaneously. 

This is not considered to detract from the technique since the cost of parts of each 
demodulator [Appendix C] is approximately $150 US [Reference 21] and can be made 
cheaper with dedicated integrated circuits. This is a significant cost reduction over an 
earlier version of the demodulator [Reference 9, page 204] which quoted the cost as $270 
US. This cost could possibly be reduced further by several means, including large quantity 
manufacture, special purpose digital signal processing (DSP) cards, or special purpose 
DSP chips. The approach of a specific DSP card follows on from the Naval Postgraduate 
School research conducted by Lt. Brian R. McGinnis, USN, and is presently being 
undertaken as a thesis project by LCDR David W. Brenner, Canadian Armed Forces. 




Figure 3.1 Multiplexing Architecture for a 3 x 4 Sensor Array. 

The possibility of reducing the number of demodulators by a system of time division 
multiplexing (TDM) of sensor information is also an option. The added complexity and 
cost of incorporating TDM may be greater than that of the extra demodulators. 
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Similarly, the relatively low cost of optical fiber for short length systems, i e , those 
where sensor, source, and receiver separation is of the order of tens of meters, has not 
made the reduction in the number of return lines a critical design parameter. There may be, 
however, an operational requirement to reduce the number of return lines. A deployable 
array or remotely operated vehicle (ROV) transmitting data to a remote location are 
possible examples where such a design consideration may be essential Methods to achieve 
significant reduction in return lines are considered later. 

An analysis of the topology of the architecture proposed in Figure 3.1 can be 
conducted to demonstrate the savings and costs associated with the multiplexing 
technique. The ratio of sensors to channels (number of optical fibers) represents a measure 
of the saving due to multiplexing that can easily be quantified. Using the architecture of 
Figure 3.1, the number of lasers is denoted by M and the number of sensors per laser by N. 
The number of fiber lines F is then given as the number of lasers plus three times the 
number of staves or return lines; 


F = M + 3A^. 


(3 1) 


Thus for the 3 X 4 array of Figure 3.1 we find the number of sensors S and number of 
fiber lines as; 


and 


S = MxN = 12, (3.2) 

F = M + W = 15, (3.3) 


S_ 

F 


MxN 
M + 3N 


0.8 sensors / channel. 


(3.4) 
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For a 10 X 10 array; 


s = 

MxN = 100, 

(3.5) 

F = 

M + 3N = 40, 

(3.6) 

100 

A/\ 

= 2. S sensors / chaimel. 

(3.7) 


Thus as the number of sensors increases the sensor-to-channel ratio increases and a 
multiplexing gain is realized. 


Table 3.1 Number of fiber lines, F, for M x A' array 


1 Sensors/laser Number of lasers M \ 

N 

2 

3 

4 

5 

6 

7 

8 

9 

10 

2 

8 

9 

10 

11 

12 

13 

14 

15 

16 

3 

11 

12 

13 

14 

15 

16 

17 

18 

19 

4 

14 

15 

16 

17 

18 

19 

20 

21 

22 

5 

17 

18 

19 

20 

21 

22 

23 

24 

25 

6 

20 

21 

22 

23 

24 

25 

26 

27 

28 

7 

23 

24 

25 

26 

27 

28 

29 

30 

31 

8 

26 

27 

28 

29 

30 

31 

32 

33 

34 

9 

29 

30 

31 

32 

33 

34 

35 

36 

37 

10 

32 

33 

34 

35 

36 

37 

38 

39 

40 


It should also be noted that this system would work with the NRL quadrature 
demodulation approach [Reference 6] which uses two of the three lines. 
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1. Power Budget Requirements. 

The limitation on the number of sensors possible is a function of the laser power 
available, the number of couplers required, and the bandwidth or frequency separation 
between modulation frequencies. The most efficient power split ratio for each coupler can 
be determined from its location in the array. The coupler ratio required for each coupler in 
a chain (horizontal row in Figure 3.1) from any laser can be calculated using Figure 3.2 as 
follows (neglecting intrinsic losses); 


A») = 

N 


(3,8) 


where P(n) is the power transmitted by coupler n to the next coupler, the output 
power of the laser, and P^ the (constant) power transmitted to each sensor The coupler 
ratio C„ is then defined as the ratio of the transmitted power, P(n), to the input power, 
P(n-l), and is given by; 


C. 


PW ^ t P. 
P(n-\) P(n-\) 


——;- for // = 0,1, ,N-\. 

+1 - // 


(3.9) 


Normalizing these calculations for a given laser output power of P„ . Equation 3.8 
becomes; 


F(//) 


P{n) 

Po 




(3-10) 


where P’(n) in Equation 3.10 represents the fraction of the total laser power available at 
coupler n. 
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Figure 3.2 Power Transmission Curve for Equal Power to Each Coupler. 

The calculations for the coupler power split ratios in of the return lines of each stave 
(vertical column of sensors in Figure 3.1) can be calculated by using the formulas from the 
previous calculations for the coupler ratios in a laser chain but for M sensors. These 
coupler ratios must then be implemented in reverse order down a stave. Thus the first or 
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top couplers in a stave of sensors will have a coupling ratio of 0.5 while the last or bottom 
couplers will have a coupler ratio of l/M. 

In a particular laser chain, power is not wasted in the sense that light is either 
coupled to a sensor or transmitted to sensors hirther along the chain. The only losses in 
the chain are splice, coupler, and fiber losses. In contrast the couplers in a stave each have 
one leg of wasted power. The most efficient use of optical power therefore will maximize 
N and minimize M The compromise is that the sensor-to-channel ratio is reduced for 
larger N and smaller M. A compromise between the two design parameters of maximum 
optical power per sensor and maximum sensor*to-channel ratio must be made based on the 
application requirements. 

A MATLAB'''^ routine [Appendix D] was written to calculate the maximum output 
power from each sensor for any given array size, M y. N. This routine takes into account 
the system losses based on the following estimates of the individual losses. These 
estimates are based on the typical losses for those components used in the 2 x 1 optical 
demonstration. Loss estimates are assumed as follows: 

0.5 dB 
0.2 dB 
3.0 dB/km 
80 m [Reference 22] 

3 mW or 4.8 dBm 
-110 dBV re 1 prad/VHz 
30 m (e g., for 50 Hz array) 


• loss/coupler (single mode) 

• loss/splice (single mode fiber) 

• fiber loss @ X = 830 nm 

• fiber length/sensor 

• typical laser output power(at Im pigtail) 

• typical noise floor at receiver 

• fiber length between sensors 

Thus the loss per sensor can be calculated as: 


2x0.5(coupler) + 4x0.2(splice) + 3x0,08km(sensorfiber) = 2.04dB/sensor, (3.11) 
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and the loss per coupler including the fiber loss between couplers as: 

0.5 + 3 X0.03km(fiber loss) = 0.59 dB. (3 12) 

The loss per coupler (0.59 dB) and loss per sensor (2.04 dB) have been included in the 
MATLAB™ routine. 

Examples are shown below for different values of M and N for both the fraction of 
laser power available and an estimated power of 3 mW, based on data for the HL8312G 
Hitachi source [Reference 23]. By comparing the calculations for a 4 x 25 element array 
and a 10 x 10 element array, the compromise between optical power and number of return 
lines for these two implementations of a 100 element array can be seen Figures 3 3 and 
3 .4 compare the fractional power transmitted to a particular coupler n in a laser chain for 
the two cases. Figures 3.5 and 3.6 compare the coupler ratios for couplers in a particular 
laser chain for the two arrays. The reason that the last coupler is not exactly one half in 
each case is due to the estimated losses included in the calculations. Figures 3.7 and 3.8 
show the power from each sensor in a laser chain for the two cases. There is 
approximately 2.7 dBm more power delivered to each sensor in the case for A^= 10. 
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Figure 3.7. Estimated power coupled to each sensor in chain for Hitachi 
HL8312G laser (3 mW at 1 m pigtail), including estimated losses, N = 25. 
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Figure 3.8. Estimated power coupled to each sensor in chain for Hitachi 
HL8312G laser (3 mW at 1 m pigtail), including loss estimation, N = 10. 





















Figures 3.9 and 3.10 show the coupler ratios for each coupler down a stave (or 
column in Figure 3.1). Again the first coupler is not exactly equal to one half as a result of 
the losses. The most important results of this analysis are shown in Figures 3 11 and 3 12 
where the estimated optical power for each sensor down a stave are compared for the two 
cases. The power for the last sensor in the stave represents the power coupled from each 
leg of each sensor in the array that is received at the photodetector It can be seen that for 
the 10 X 10 array the power lost as a result of the lost (unused) power from each coupler 
in every stave is so significant that there wilt insufficient power at the photodetector. The 
estimated system noise floor is expected to be approximately -90 to -110 dBm The 4 x 
25 array, on the other hand, is able to deliver approximately -58 dBm, which is on the 
order of 50 dB above that for the 10x10 array. 
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Coupling Ratio 
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Coupler #(M.1 total) 

Figure 3.9. Coupler ratio for couplers in each stave (column in 
Figure 3.1) for 4 x 25 element array. 
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Figure 3.10. Coupler ratio for couplers in each stave (column in 
Figure 3.1) for 10 x 10 element array. 






















An optical sensor array topology is not only dilTerent from present piezoelectric with 
twisted pair, or coaxial lines, but may vary depending on the specific application. A hull 
mounted array, for example, may best be implemented with no multiplexing in either one 
laser per sensor, or a shared laser topology as shown in Figure 4.1. The savings in size and 
weight realized by an optical system over those presently available may preclude the 
requirement for multiplexing 




Figure 3.13. Array topologies for a hull mounted optical array: 
one sensor per laser (upper) and shared laser chain (lower). 
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2. Bandwidth Limitations 

Consider a system of M lasers with the lowest modulation frequency j at 20 
MHz. Let the upper design limit of 20 kHz be the maximum acoustic frequency with a 
dynamic range of 140 dB at 100 Hz Using Equation 2.10, and assuming similar noise 
floor and maximum acceptable distortion levels, this gives the significant bandwidth of the 
interferometric signal as; 

A = = 5x(0.5x10^x10““'“)x 20 = 500 kHz (3.13) 

With the fact that we have both upper and lower sidebands, a 1 MHz frequency separation 
is required between laser modulation frequencies. For at 20 MHz this represents an 
upper limit of 20 sensors that can be accommodated before the second harmonic of the 
lowest carrier . 

Investigation of an M x array reveals that the optical power restrictions imposed 
by coupling sensors in a stave and the subsequent lost or wasted power at every coupler is 
a far more stringent design consideration than bandwidth limitations 

B. REDUCTION IN THE NUMBER OF RETURN LINES. 

An application such as a deployable array or remotely operated vehicle (ROV) may 
require significant reduction in the number of return lines. The number of return lines may 
be required to be reduced beyond that shown in Figure 3.1 or even down to a single return 
line. The cost of additional optical fiber may be greater than that of the extra couplers and 
processing required in a system with reduced return lines. This is the case in systems 
where the distance between sensors and/or receivers and/or sources become too large. At 
approximately 50 mik iv the gallon (container), fiber is far cheaper and more co'np.'ct 
than copper. 
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A simple method for reducing the number of return lines utilizing 3x3-coupler 
terminated sensors and symmetric demodulation is a combination of OSCUM detailed 
earlier and time division multiplexing (TDM). TDM has been successfully demonstrated 
with intensity based sensors [e g.. Reference 2, pages 165-6, and Reference 4, page 95] 
and interferometric sensors [e g.. Reference 2, page 172-176, Reference 3, pages 82-84, 
and Reference 4, pages 96-100]. 

The discussion of Reference 3, page 83, in which TDM with Mach-Zehnder sensors 
terminated with 2x2-couplers is presented, can be adapted to the case for 3x3-coupler 
terminated sensors and symmetric demodulation. Such a scheme is shown in Figure 3.14. 
The demultiplexing of the pulses would have to be carried out prior to demodulation 
Each input pulse generates one pulse per sensor on each of the three output lines. These 
output pulses must be appropriately gated to relate them to their respective sensor In 
order to correctly determine which pulse is associated with which sensor and to ensure 
that there is no overlap of pulses, the pulse duration must be less than the propagation 
delay between any two sensors, and the time between pulses must be greater than the 
propagation delay difference between the shortest and longest paths. 

The power output per sensor can be maximized in a fashion similar to the previous 
section by choosing the coupling ratios of the couplers appropriately. There is, however, 
an additional power loss in a TDM system due to the duty cycle of the source. 
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Figure 3.14. TDM willi 3x3-cou|)ler terminated sensors For symmetric 
demodulation (adapted from Reference 3). 


1'he number of return lines of the complete array could be reduced to a single line by 
1 DM of the individual sensor legs. This introduces further losses and demultiplexing 
complexity. A combination of OSCUM and TDM could be employed to improve the 
optical power per sensor (equivalently the SNR of each sensor). A scheme incorporating a 
single return line and a combination of OSCUM and TDM is proposed in Figure 3.15. The 
loop diameters indicate relative differences in path length delays. 
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rigiire 3.15. Single rcluni arcliitcclure iiicorporaling OFSUM and I'DM. 

It sliould also be noted that a system with a reduced number of return lines is not as 
robust in the sense that the loss of a single fiber or coupler leads to a much larger or total 
loss of sensor information while this loss is reduced for the case of more return lines. 
Again the application will dictate the costs and benefits involved. References 3, 4. and 5 


57 










provide an overview of the techniques presently used with a variety of other types of 
sensors and demodulation techniques. Further research is required to apply these 
techniques to 3x3-coupler terminated Mach-Zehnder sensors with symmetric 
demodulation. 
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rv. CONCLUSIONS 


The technique of 'optical fiber subcarrier frequency multiplexing' has been successfully 
demonstrated with a 2 x 1 element array. This is the first time that this multiplexing 
technique has been known to be demonstrated with interferometric sensors. 

The architecture for a particular array is very application specific with the requirements 
of the particular application dictating the multiplexing technique to be utilized. This occurs 
as a result of the large power budget trade-offs incurred when optical signals are 
recombined through couplers. The development and introduction of low cost optical 
amplifiers may provide a solution to this problem. 

The following list details some of the follow-on work recommended to further the goal 
of all optical sensor systems for a wide variety of applications. 

• Demonstration of the technique at RF frequencies with a larger number of 
sensors (or simulated sensors), possibly with highly linear sot a as similar to 
those described in Chapter II, section B.2. 

• Investigation of a method to digitally implement the multiplexing scheme to 
further reduce costs and hardware complexity. 

• Investigation and demonstration of techniques to reduce the number of 
return lines for specific applications. 

• A less expensive 3x3-coupler, or alternative to the 3x3-coupler that 
presently terminates the sensors developed at NFS. 
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APPENDIX A. SPECTRUM OF AN INTERFEROMETRIC SIGNAL 


The interferometric signal output from a 3x3 terminated Mach-Zehnder interferometer 
with input signal <f)(t) (Equation Al), expanded into a Bessel weighted series, is given by 
equation A2; 


I(t) = A + Bcos(^t)) 
^t) = ^.sin<u,t. 


(Al) 




Jo(^s) + 2Y,JM)cosi2lco^t) 


/=! 




21 -/ 2/»l (^,)sin [( 2 / + l)<y,/] 


»=o 


(A2) 


where ^jri/i(0 is the slowly varying (compared to the signal of interest) or drift phase term 
which includes the phase shift due to environmental effects (hydrostatic pressure, 
temperature, etc.) as well as the physical difference in the length of the legs and the phase 
shift due to the coupling region of the interferometer. 

Thus, for a given (assumed to be a constant over the time scale of the signal 

of interest 27t/<y^) the ratio of the spectral components of odd or even harmonics of the 
signal of interest will be determined by the ratio of the Bessel functions. This property can 
be used to determine the operating point, i.e., the magnitude of by measuring these 
ratios of components and relating them to the equivalent Bessel function ratios (i.e., J 2 /J 4 , 
J 3 /J 5 , etc.) [Reference 9, page 28]. This is especially useful for determining the operating 
point in the 'sub-fnnge' range. Figure AI shows the lower order Bessel functions and 
Table Al, the component ratios for = n radians amplitude, which can be used to verify 
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the examples of * Ji radians amplitude and drift angles ^dnfi ~ 120®, and 240®, 

shown in Figure A2. 













































































APPENDIX B. SYMMETRIC ANALOG DEMODULATOR VERSION 2.1 


Reference 9 details the first implementation of the symmetric demodulation scheme 
using analog electronics. The analog implementation was redesigned by Dr. David L 
Gardner with vastly improves the small-signal performance over the original 
implementation. Gardner's implementation was not previously documented and was used 
for the multiplexing demonstrations of this thesis. 

The circuit schematic is shown at Figure Bl. The algorithm is well documented 
[References 3, 4, 5, and 7] with version 2.1 described here being configured for lower 
noise floor and ease of component manipulation. The circuit description is divided into 
several distinct stages adapted from Dr. Gardner's notes. Performance results of version 
2.1 are included. 

A. INPUT STAGE 
1. Photodiodes. 

Three input choices are available as follows. 

1. Three 1300 nm photodiodes are the lower row of bulkhead SMA connectors 
Presently fitted are model UDT080 devices (with an active area of 0.08 mm^). These 
devices are selected by opening the box and moving each of the three blue jumpers 
(adjacent to the edge connector) to the left-most position as one looks down on the board 
facing the edge connector. These small and, hence, fast diodes may present some difficulty 
in alignment. Both the fiber and diodes can be epoxied into position for field tests. 
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2. Three 830 nin photodiodes are the upper row of bulkhead SMA connectors 
Presently fitted are model HRS 102 devices (with an active area diameter of 300 pm) 
These detectors are selected by opening the box and moving each of the three blue 
jumpers (adjacent to the edge coimector) to the right-most position as one looks down on 
the board facing the edge connector. These also should be fixed into position with epoxy 
for field tests. 

3. It is often convenient to use interferometric sensor simulators to test 
demodulators and conduct various other trials. Since the output fi-om these simulators is a 
voltage, an input resistor has been installed to allow the demodulator to be used with the 
simulators. This will cause the input amplifier to be an inverting type rather than the 
transimpedance configuration required for the photodiodes Three 9207 Q resistors are 
precision type and can be used in conjunction with the fixed value feedback resistor to test 
the operation of the demodulator input. The box must be opened to gain access to these 
resistors. They are gray and installed vertically along the edge of the board adjacent to the 
connector. The long lead of the resistor is available for connection but, first, the blue 
jumpers must be removed. 

2. Input Amplifiers. 

Two options are available for gain control of the input amplifiers: 

1. Fixed-value resistors - In the course of testing interferometric sensors, it is 
useful to determine the light power in each of the sensor fibers. The typical responsivity of 
the photodiodes is 0.5 AAV. Since the current-to-voltage ratio or conversion factor must 
be known to make such a measurement, 20 kH resistors have been installed. These 
resistors are selected by opening the box and moving the black jumpers, located along the 
left edge of the board as viewed facing the edge connector, to the leftmost position. 
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2. Potentiometer-controlled gain - The output power from each leg of a coupler 
is rarely uniform. Variations in the modulation depth, splice losses, bending losses and 
modulation losses all reduce the light emitted from each fiber. Since the demodulator 
algorithm assumes equal power from each leg, potentiometers have been installed to 
provide adjustment of the amplitude of the AC portion of each input signal. Adjustments 
of the DC portion may also be accomplished. The potentiometers are blue in color and are 
located just in from the left-hand edge of the board as viewed facing the edge connector. 
The output of each amplifier may be monitored by connecting an oscilloscope lead to the 
test point loops near the yellow jumpers and aligned with the transverse axis of the board. 
These test points are connected to pin 5 of the non-inverting amplifiers (ICs U4-U6 of 
Figure Bl). 

B. INTERMEDIATE AMPLIFIERS AND OUTPUT BUFFERS 

1. Intermediate Amplifiers 

The intermediate amplifiers have been included to permit adjustment of the signal 
level to the remainder of the circuit. The gains of these non-inverting amplifiers are 
selected via the yellow jumpers towards the left edge of the board, as viewed facing the 
edge connector. The gains permitted are 20, 30, and 40 dB, using pairs of jumper leads, 
right to left as viewed facing the edge connector. The signal levels at this point in the 
circuit may be monitored using the test loops adjacent to these intermediate amplifiers, 
aligned with the longitudinal axis of the board. 

2. Output Buffer Amplifiers 

For convenience in recording and monitoring requirements, the intermediate 
amplifier outputs are also connected via buffer amplifiers (U8A, USB, and U8C in Figure 
Bl) to the lower row of BNCs available on the front panel. These correspond to channel 
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numbers 1-3 (as numbered on the left hand edge of the board) from left to right, viewed 
facing the front panel. 

C. MODULATION DEPTH ADJUSTMENT 

The modulation depth detected in each fiber has been observed to vary slightly 
between legs. Adjustment of the AC portion for equal signal level between channels also 
changes the DC portion of the signal. A set of summing amplifiers has been included to 
remove the DC portion of the signal. This approach an alternative to high-pass filtering the 
input signal since sub-radian interferometric signals contain DC. This approach has proven 
successful, provided that the ratio of modulation depth variations between channels 
remains constant. Changes in the ratios of 3x3-coupler outputs seems to cause problems 
Investigation of the demodulator performance is required to address interchannel 
variations. 

The adjustment to remove the DC portion of the signal is made via the potentiometers 
(tan colored) near the middle of the board. The effect of this adjustment to each channel 
can be monitored using the test loops under the differentiator resistors (mounted on 
component towers in the middle of the board). The signals at these test points should be 
adjusted such that they have equal AC amplitudes and zero DC component. Adjustment to 
correct discrepancies should be performed using the input amplifier potentiometers (blue) 
first and then the modulation depth potentiometers (tan). 

D. DIFFERENTIATOR ADJUSTMENT 

The signal presented to the multiplier sections should be as large as possible without 
distortion (10 V). To this end the differentiator transfer function at the highest frequency 
of interest should be unity. The resistors mounted on the component headers permit 
adjustment of this transfer function. In this implementation 40.2 kO resistors are fitted. 


69 






corresponding to an upper frequency of interest of approximately 40 kHz, as computed in 
the following formula. 


1 

2^RC 


_ 1 _ 

2;r(40.2knX100pF) 


39 6 kHz, 


(Bl) 


E. PERFORMANCE 


Reference 9, Chapter XI, details the performance of the first implementation of the 
symmetric analog demodulator (SAD); Reference 9 also details a digital implementation of 
an Asymmetric Digital Demodulator. Following the analysis of Reference 9, page 161, the 
analog interferometric simulators were used to provide inputs to the SAD V2.1. This 
removes the added noise due to laser phase and the wandering or drifi of the polarization 
of the optical waves in the interferometer, the major cause of drift in amplitude of the 
interferometric signal. The Scale Factor of the demodulator is defined as the ratio of the 
demodulator output voltage modulation amplitude to the phase modulation amplitude of 
the input interferometric signals. 

The phase noise for the demodulator at 5 kHz was measured as 6.33 prad/VHz The 
average phase noise spectral density of the three analog interferometric simulators was 
measured at 1 kHz to be 0.48±0.005 prad/>/Hz [Reference 15, book 1, page 19] and 
although not measured at 5 kHz is still far enough below the phase noise floor of the 
demodulator that it would not be of the order of the demodulator noise at 5 kHz. Thus the 
demodulator noise measurement can be mostly attributed to the demodulator itself The 
dynamic range at 5 kHz was determined for maximum signal (4% THD) to be; 

= 7.4x10* 

6.33/a-ad (B2) 

Dynamic Range at 5 kHz = 201og(7.4x 10*) = 137 dB. 
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APPENDIX C. DEMONSTRATION WITH ANALOG ELECTRONICS 


Prior to an optical system demonstration a single channel demonstration of the 
multiplexing technique was attempted for proof of concept purposes. Since a 
demonstration with only a single interferometric output from the 3x3-coupler was 
attempted, the demodulation of the interferometric signals could not be carried out. The 
essence of this experiment was that if the ratios of the demultiplexed components of the 
interferometric signal and the original signal remained constant throughout the various 
steps of frequency translation, then one could conclude that the technique was feasible 
That is to say, if the fundamental and significant higher order harmonics, which are greater 
than the noise floor, were in the same ratios as the original interferometric signal, 
demodulation would produce the same output phase modulation for both the baseband 
and the demultiplexed signals. Further, the demonstration did not quantify the effects of 
the loss of higher order harmonics (as a result of the multiplexing scheme). The effect 
would be extremely small but could be quantified in terms of the amount of energy lost 
between the baseband signal and the demultiplexed signal. This is easily computed from 
the fact that the energy contained in the interferometric signal that can be demodulated is 
the sum of the amplitude squared values of all of the components (Parseval's theorem). By 
calculating the energy before and after multiplexing, and the energy lost, a reasonable 
estimate of the effect of the distortion fi'om the demodulated signal can be made. It will be 
shown later that this effect is small. For low level signals however, it provides a lower 
limit on the level of signal which can be demodulated, i.e., the minimum detectable signal. 
This analysis assumes that dispersion of the higher order harmonics is minimal, i.e., the 
higher order harmonics are integral multiples of the fundamental. 
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Another aspect of the demultiplexing technique not covered in this demonstration is 
the effect of filtering on the interferometric signal prior to demodulation This was not 
considered since filter theory is well understood and documented, and the only 
requirements to ensure successful demodulation are constant gain and phase over the pass 
band of interest. This passband of interest is the interferometric signal significant 
bandwidth/,^ discussed earlier (section II B I .). 

The equipment setup for the demonstration is shown at Figure Cl. 
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Figure Cl. Equipment Setup for Aiinlog Electronics Dciuoiistrntion. 
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The setup has three separate function generators to simulate the acoustic signal from 
the sensor and wavelength and intensity modulation from the laser source The wavelength 
modulation is phase and frequency locked to the intensity modulation of the source via 
appropriate operation of the triggers of the respective function generators. Another 
function generator is utilized as the demultiplexer or mixer by appropriate use the AM 
input. The outputs were monitored by a HP3562A spectrum analyzer. 

A. ANALOG INTERFEROMETRIC SIMULATOR 

The analog interferometric simulator is similar to that described at Reference 9, 
Appendix C. A modification was conducted so that the input previously labeled SENSOR 
2 is now labeled STATIC PHASE 0/P and provides an easier way to measure the static 
phase drift angle Previously this was determined using Lissajous figures [Reference 
20, pages 30-39] and making appropriate measurements [Reference 9, Appendix C, 
section B, pages 301-2]. The drift angle term in the simulator is not time-varying as it 
would be in an actual system but a constant set on the front panel via the static phase 
adjustment. 

The modified simulator block diagram is shown in Figure C2 while the circuit is 
shown in Figure C3. The static phase measurement calculations will be presented by 
isolating that part of the circuit involved, see Figure C4 
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Figure C3. Modified Analog Iniviferometric Signal Simulator Circuit Diagram. 
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Tlie Analog Devices AD639 Universal Trigonometric Function Converter is 
conllgurcd to produce the sine of its input. A one-volt input is equivalent to an input 
argument of 50° Thus the gain G of the circuit of Figure C4 is: 


G 


2 


2.85. 


(Cl) 


Thus the static phase adjustment input equates to; 

2.85 x 50°/Volt = l42°/Volt. (C.2) 


so that, with one simulator adjusted to zero volts static phase output and the other two 
static phase outputs adjusted to ±0.845 V, we obtain our three interferometric signals 
symmetrically phase shifted 120° relative to each other. 


STATIC PHASE 
MEASUREMENT 



Figure C4. Circuit for Static Phase Measurement Calcuiatioiis. 
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B. RESULTS 

The first test verified that the modulated components were in the same ratio as the 
baseband components. The input labeled PZTMOD in Figures Cl and C2 was grounded 
to simulate the absence of wavelength modulation. The static phase was set at 120° 
(0.845V static phase output) to ensure that both even and odd harmonics would be 
present in the interferometric signal. A 1 kHz sine wave was applied to the SENSOR 1 
input to generate an interferometric type signal The amplitude was varied and the ratios 
of baseband components and modulated components were observed to 'track' each other 
A plot of a4.2Vpp amplitude = 2.1 radians from simulator scale factor of 1 rad/V) is 
shown at Figure C5 with the upper waveform showing the baseband 'interferometric 
signal' prior to multiplexing, and the lower signal showing the amplitude modulated and 
demultiplexed version. 
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Figme C5. Speciruiti of baseband simulated interferometric signal (top) and 
speciriiiii of demultiplexed version of same interferometric signal 
(bottom). 







It can be seen that signal and demultiplexed components are indeed in the same ratios 
even though the amplitudes are different. 

Figure C6 shows the Lissajous pattern for the demultiplexed output (x axis) and the 
simulated interferometric signal output (y axis) with no wavelength modulation (PZT 
MOD input grounded). This also verifies that the signals are in the same phase and of the 
same shape, i.e , the component ratios are the same. The closed ellipse or nearly straight 
line offset at 45° indicates a zero phase difference. The total harmonic distortion THD is 
the same for both waveforms. The THD together with Figure C5 also provides a crude 
measure of the ratios of components. 

The HP3562A Spectrum Analyzer calculates the THD by calculating the ratio of 
energy in the fundamental component to the sum of the other components Thus while it is 
possible to have the components in different ratios and the same THD, it is not possible to 
have them in different ratios and the same THD. 
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Wavelength modulation of the source is simulated by a variable amplitude signal 
locked in phase and frequency to the modulation frequency connected to the PZT MOD 
input. This wavelength modulation was not included in Equation 13 (where it was 
presumed negligible). Equation 13 is rewritten below with the wavelength modulation 
included. 


u = A + Bcos(^t)+^j^(t)+^,+^„sinn>„t). (C.3) 


Hence represents the amount of wavelength modulation. Figure C7 shows the 
Lissajous pattern for the case where (j>,„ = 0.42 radians for an 'acoustic' signal of 2 12 
radians amplitude at 1 kHz signal and a modulation frequency yj,, = 20 kHz This clearly 
shows that the two signals are no longer in phase (open ellipse) or the same shape (seen as 
the jagged shape of the outside of the ellipse). This is further highlighted by Figure C8 
where the THD ratio and component ratio are no longer the same 

This change highlights the difference between this multiplexing scheme and the 
presently popular wavelength modulated frequency division multiplexing (WM-FDM) 
WM-FDM is more commonly referred to as phase generated carrier (PGC) [Reference 5, 
section 11.4.2.2, pages 360-2 and Reference 1, volume II, section 14.6.5, pages 564-5]. 
While WM-FDM seeks to maximize the wavelength modulation and to reduce the 
intensity modulation, our technique AM-FDM operates optimally with no wavelength 
modulation. 

The analog electronics demonstration showed that the multiplexing technique will 
work provided that wavelength modulation is not significant. The extent to which the 
proposed multiplexing technique will work with increasing or decreasing wavelength 
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Figure C8. Spectriiiii of signals of Figure C7 (with wavelength modulation); 

baseband interreronietric signal (top) and spectrum of demultiplexed 
version of same inteiTerometric signal (bottom). 
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APPENDIX D. MISCELLANEOUS 


This appendix contains the following; 

• MATLAB™ routines written by the author for this thesis. 

• THB292 22 kHz integrated active filter response. 

• Operating curves and data sheets for lasers used in this thesis. 

• Benchtop Interferometer Figure. 

• ORTEL™ information sheet on device linearity specifications. 

• Noisefloor Figures for Analog Interferometric Simulator and Mixing circuit 

The MATLAB™ routines are provided with no guarantee of results obtained i e., the 
author disclaims any errors in the code reproduced. 
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MATLAB™ routine used for Chapter II figures. 


% Thciis - Mixing of intoferomcthc signal simulation 
clear.clg; 

61»lOOO^ihis 1 

62-600ailus2”pi;phkl2-0; 
finl>40000; 
fin3>60000; 
phis>pi; 

% phid is the drift phase term for the difference in optical path length 

% between the two legs of the uiterferomeler 
%iphid>input('£nter phid (in radians) *); 
t-0:( 10/(61 •32767));(10/&1); 

%l=0;( 10/(fs»32767));( 10/fs); 
m=0.75. 

11«(1 +(cos(phisl •sin(2*pi*61 •t) +phidl ))).*( 1 +(m«sin(2»pi»lmi»t))); 
I2=( 1 +(cos(phis2*sin(2*pi*fs2»t) +phid2))).*( 1 ^^m*sin(2*pi»fin2»l))); 

I3=(n-<-I2); 

I=I3.*sin(2*pi*fml*t); 

F=ffl(1.32768). 

f=0; l/iiiax(t):32767* l/max(t). 

(k=f./1000; 

subplot(21 IXplotIff (1:1200),20*logl 0(abs(F( 1; 1200)))).grid 
xlabeK'FrequeiKy (kH 2 )'),ylabelCMagnitude (dB)') 
axis((0 fled 200)0 100)); 
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MATLAB™ routine used for Chapter HI figures. 


SThcai* nlculatiom for MxN array 
clear,ctg; 

N-inpitfCEiiter # of icnaon N *); 
M'lnputCEnlcr * of icmaci M 7i 


fori-l:N 


I(iH( I -1 • -(i)/N>+ I/Ny J .230268770812; 

e(i)=(( 1 - l/NyNy(( 1 - 1/N)»( 1l/N); 
pti(i)=l(i)*c(i); 

pkl(iH).003*I(i)*e(i); 

end 

c' 

fdB= 10*logl 0(anes(N, 1 )-c’) 

pkdb= 10*log 10(pk);pk 1 db*’10*log 10(pk 1), 

%c I ^flipud(c< I :lengUt(c>-l)'). 

<^JI(l)-pk(l); 

•/«forj=2:lenglh(c>-l 

O/Wnd 


for i=l:M 


InKiHC 1 -1 /M)*( KiyM)-^ 1/My 1 ■ 145512941446; 

cni(i)=((l-l/MyMy((l-l/M)*(l-i/M)*l/M); 

•/^(i)=l(i)*c(i); 

•/ipkl(i)=0.003»l(i)*c(i); 

end 

cin' 

fdBin= 10*logl0(ones<M, 1 )<ni') 

%pkdb= 10*iogl0(pk);pk I db= 10’logl 0(pk I); 
c 1 =f]ipud(cn)( 1 :length(cni)-1)'); 
ll(l)=plcl(l); 


for i-2:lenglh(cm) 

Il(i)-(cl{i-l)Ml(i.l)); 

end 

IIdB°10*logl0(Il); 

figurefl) 

plol(I,'o7 

ylabel(Tractional Power Input Pin/Ptotal') 
vlabelCSentor # N'Xgrid 


figuirv(2) 

plo«(c,V).yUbeK‘Coupl«' RaJio) 

KlabeK'Coupkr »(N-1 toul)').gn<i 

sxu((ON 0 0.61) 

fif{ure<3) 

plo<(fda'o'),yUbel(X^pta' Raiio (dB)') 
xIabeU'Coupler tt (N-1 U)Ul)’).grid 
axis((0 N -5 01) 

figure(4) 

plot(pkdb,'o‘) 

yUbel('Pow«r to Sensor N Pn/PloUl (dB)') 
xIabelCScnsor *» N‘),grid 

figure(}) 

plol(pkldb-2.04,'o'),ylabel('Pow«f Sensor N Po“3niW (dBm)') 
xlabeK'Sensor 4 N'),grid 

rigure<6) 

plol(cl.'o'),grid 

ylabeICCoupliiig Ralio'),xlabe)CCoupler » (M-1 toul)') 
figure(7) 

plol(l I dB-2.04,'o‘),gnd 

ylabel('Power at Sensor M (dB) ),xlabel('Sensor« M') 



MATLAB™ routine used for figure A2. 


% IiHetfcronielhc Signal Simulation 
S 

H for A/B= 1 1(1)-1 +ooi( phis (sin(oniegai l] + phid) J 
% 

clcar.clg; 

&-I000; 

% phis is the amplitude of the signal of interest 

phis-input(^ter phis (in radians) *); 
Sphisl-inpulCEnter phisl (in radians) *); 
Sphis2-input('Enter phis2 (in radians)'); 


% phid is the drift phase term for the difference in optical path length 
% between tlie two legs of the interferometer 
phid-input('£nter pliid (in radians) ‘); 
phidl-input('Enter phid (in radians) '); 
phid2-input(‘Enter phid (in radians) ‘); 

t=0.(l0/(fs*32767)).(10/fs); 

1= I ■t-(cos(pliis*sin(2*pi*fs*t) -rphid)); 

11 = 1 +(cos(plus*sin(2*pi*fs»t) -t-pltidl)), 

12=1 +(cos(pliis*sin(2*pi*fs*t) •>-phid2)); 

F=fll(1,16384); Fl=ffl(11.16384); F2=ffl(12.16384); 
f-0:2/max(t); 16383*1 /max(l); 

subplot(221),plot(t(l:length(t)/1000:lengtl*'t)'5),l(l:lenglh(t)/1000:length(t)'5)),grid 
xlabel(Time (sec)'),ylabel('lntensit>’‘) 
axis([01(70)0 100|); 

subplot(222),plot(l(1.400),20*logl0(abs(F(l:400)))).grid 
xlabel(Trequency (Hz)'),ylabel(‘Magnitude (dB)') 

$ubpIot(223),plot(t(I:length(t)/l000:lengtli(t)/5),II(l:length(t)/1000:length(t)'S)),grid 
xlabel(Time (sec)'),ylabel(’lnlensity’) 
axis((01(70)0 100]); 

subplot(224).plot(f( 1:200),20*logl 0(abs(F 1 (1:200)))),grid 
xlabelCFrequency (Hz)'Xylabel('Mapiitude (dB)') 
pause 
cig; 

subplot(22I),plot(t(l:lenglli(tyiOOO:length(ty)),I2(t:length(iyiOOO:length(t)/5)),grid 

xlabelfTime (secyXylabelCIntensity*) 

pause 

axis([0 f(70)0 100]); 

subplot(222Xplot(i( 1:4000),20*logl 0(abs(F2( 1:4000)))Xgrid 
xlabelCFrequency (Hz)'),ylabelCMagnitude (dB)') 
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Figure D2. Seaslar Laser #SP2676 Maiiiiractiirer's lest data. 
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Figure D3. Seastar Laser #SP2676 measured operating curve. 
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Figure D5. Seaslar Laser //SP1320 measured operating ciii've. 
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Figure D7, ORTEL™ device linearity speciricalions. 
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